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Abstract — Tracking Lagrangian coherent structures in dy- 
namical systems is important for many applications such as 
oceanography and weather prediction. In this paper, we present 
a collaborative robotic control strategy designed to track stable 
and unstable manifolds. The technique does not require global 
information about the fluid dynamics, and is based on local 
sensing, prediction, and correction. The collaborative control 
strategy is implemented on a team of three robots to track 
coherent structures and manifolds on static flows as well as a 
noisy time-dependent model of a wind-driven double-gyre often 
seen in the ocean. We present simulation and experimental 
results and discuss theoretical guarantees of the collaborative 
tracking strategy. 

I. INTRODUCTION 

In this paper, we present a collaborative control strategy 
for a class of autonomous underwater vehicles (AUVs) to 
track the coherent structures and manifolds on flows. In 
realistic ocean flows, these time-dependent coherent struc- 
tures, or Lagrangian coherent structures (LCS), are similar 
to separatrices that divide the flow into dynamically distinct 
regions. LCS are extensions of stable and unstable manifolds 
to general time-dependent flows [1], and they carry a great 
deal of global information about the dynamics of the flows. 
For two-dimensional (2D) flows, LCS are analogous to ridges 
defined by local maximum instability, and quantified by local 
measures of Finite-Time Lyapunov Exponents (FTLE) [2]. 

Recently, LCS have been shown to coincide with optimal 
trajectories in the ocean which minimize the energy and the 
time needed to traverse from one point to another [3], [4]. 
Furthermore, to improve weather and climate forecasting, 
and to better understand various physical, chemical, and 
geophysical processes in the ocean, there has been significant 
interest in the deployment of autonomous sensors to measure 
a variety of quantities of interest. One drawback to operating 
sensors in time-dependent and stochastic environments like 
the ocean is that the sensors will tend to escape from their 
monitoring region of interest. Since the LCS are inherently 
unstable and denote regions of the flow where more es- 
cape events may occur [5], knowledge of the LCS are of 
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paramount importance in maintaining a sensor in a particular 
monitoring region. 

Existing work in cooperative boundary tracking for robotic 
teams that relies on one-dimensional (ID) parameteriza- 
tions include [6], [7] and [8], [9] for static and time- 
dependent cases respectively. Formation control strategies for 
distributed estimation of level surfaces and scalar fields in 
the ocean are presented in [10]— [12] and pattern formation 
for surveillance and monitoring by robot teams is discussed 
in [13]— [15]. Our work is distinguished from existing work 
in that we use cooperative robots to find coherent structures 
without requiring a global picture of the ocean dynamics. We 
take inspiration from [16] and design a strategy to enable 
a team of robots to track the stable/unstable manifolds of 
general 2D conservative flows through local sensing alone. 
We verify the feasibility of our method through simulations 
and experiments and show how the proposed strategy can 
be extended to track coherent structures in time-dependent 
conservative flows with measurement noise. To our knowl- 
edge, this is the first attempt in the development of tracking 
strategies for mapping LCS in the ocean using AUVs. 

The novelty of this work lies in the use of nonlinear 
dynamical and chaotic system analysis techniques to derive 
a tracking strategy for a team of robots. The cooperative 
control strategy leverages the spatio-temporal sensing capa- 
bilities of a team of networked robots to track the boundaries 
separating the regions in phase space that support distinct 
dynamical behavior. Additionally, our boundary tracking 
relies solely on local measurements of the velocity field. 
Our technique is quite general, and may be applied to any 
conservative flow. 

The paper is structured as follows: We formulate the 
problem and outline key assumptions in Section III] The 
cooperative control strategy is presented in Section |IlT] and 



its theoretical properties analyzed in Section IV Section 
|V| presents our simulation and experimental results. The 
extension of the proposed strategy to a noisy time-dependent 



model of a wind-driven double-gyre is presented in VI We 



conclude with a discussion of our results and directions for 



future work in Sections VII and VIII respectively. 



II. PROBLEM FORMULATION 

We consider the problem of controlling a team of N 
planar AUVs to collaboratively track the material lines that 
separate regions of flow with distinct fluid dynamics. This is 
similar to the problem of tracking the stable (and unstable) 
manifolds of a general nonlinear dynamical system where 
the manifolds separate regions in phase space with distinct 
dynamical behaviors. We assume the following 2D kinematic 



model for each of the AUVs: 

Xi = Vi cos 6j + iii, (la) 

y i = V i smd i +v i , (lb) 

where x,- = [xj, y t ] T is the vehicle's planar position, V, and 0, 
are the vehicle's linear speed and heading, and u,- = [«,-, vi\ 
is the velocity of the fluid current experienced/measured by 
the i th vehicle. Additionally, we assume each agent can be 
circumscribed by a circle of radius r, i.e., each vehicle can 
be equivalently described as a disk of radius r. 

In this work, u,- is provided by a 2D planar conservative 
vector field described by a differential equation of the form 

x = F(x). (2) 

In essence, m, = F r (x,) and v, = Fy(xj). Let B$ and By denote 
the stable and unstable manifolds of p}. In general, B$ and 
Bu are the separating boundaries between regions in phase 
space with distinct dynamics. For 2D flows, B* are simply 
one-dimensional curves where * denotes either stable (S) or 
unstable (U) boundaries. For a small region centered about a 
point on B*, the system is unstable in one dimension. Finally, 
let p(B*) denote the radius of curvature of B* and assume 
that the minimum of the radius of curvature p, m „(B*) > r. 
This last assumption is needed to ensure the robots do not 
lose track of the B* due to sharp turns. 

The objective is to develop a collaborative strategy to 
enable a team of robots to track B* in general 2D planar 
conservative flow fields through local sampling of the veloc- 
ity field. While the focus is on the development of a tracking 
strategy for B$, the proposed method can be easily extended 
to track Bj/ since Bu are simply stable manifolds of |2]i for 
t < 0. We present our methodology in the following section. 

III. METHODOLOGY 

Our methodology is inspired by the Proper Interior Max- 
imum (PIM) Triple Procedure [16] - a numerical technique 
designed to find stationary trajectories in chaotic regions with 
no attractors. While the original procedure was developed for 
chaotic dynamical systems, the approach can be employed to 
reveal the stable set of a saddle point of a general nonlinear 
dynamical system. The procedure consists of iteratively 
finding an appropriate PIM Triple on a saddle straddling 
line segment and propagating the triple forward in time. We 
briefly summarize the procedure in the following section and 
refer the interested reader to [16] for further details. 

A. The PIM Triple Procedure 

Given the dynamical system described by |2]), let ®€l 2 
be a closed and bounded set such that 2> does not contain 
any attractors of |2]). Given a point xe^, the escape time of 
x, denoted by 7e(x), is the time x takes to leave the region 
S> under the differential map Q. 

Let J be a line segment that crosses the stable set B^ 
in 3>, i.e., the endpoints of the J are on opposite sides of 
B$. Let {x£,xc,Xfl} denote a set of three points in J such 
that xc denotes the interior point. Then {xl,xc,xr} is an 



Interior Maximum triple if 7e(xc) > max{7£(x£), Te(xr)}. 
Furthermore, {xi,xc,xr} is a Proper Interior Maximum 
(PIM) triple if it is an Interior Maximum triple and the 
interval [xl,xr] in J is a proper subset of J. 

Then the numerical computation of any PIM triple can 
be obtained iteratively starting with an initial saddle straddle 
line segment Jq. Let x^ and x# denote the endpoints of Jq 
and apply an £0 > discretization of Jq such that x^ — qo < 
qi < ... < qM = x/; . For every point q„ determine 7i(q,) 
by propagating q, forward in time using Q. Then the PIM 
triple in Jq is given by the the points {q,t_i , q^, q^+i } where 
qk = arg max 7k (q,). This PIM triple can then be further 

i=l,.. .,M 

refined by choosing J\ to be the line segment containing 
{Qjt-ijQ/tjQjt+i} an d reapplying the procedure with another 
£1 > discretization where £1 < £0. 

Given an initial saddle straddling line segment Jq, it has 
been shown that the line segment given by any subsequent 
PIM triple on Jq is also a saddle straddling line segment [16]. 
Furthermore, if we use a PIM triple x(f) = {xl,xc,xr} as 
the initial conditions for the dynamical system given by |2]) 
and propagate the system forward in time by At, then the 
line segment containing the set x(f +Af), Jt+At, remains a 
saddle straddle line segment. As such, the same numerical 
procedure can be employed to determine an appropriate PIM 
trip on J t +At- This procedure can be repeated to eventually 
reveal the entire stable set B$ and unstable set By within S 
if time was propagated forwards and backwards respectively. 
Furthermore, since the procedure always begins with a 
valid saddle straddling line segment, by construction, the 
procedure always results in a non-empty set. 

Inspired by the PIM Triple Procedure, we propose a 
cooperative saddle straddle control strategy for a team of 
N > 3 robots to track the stable (and unstable) manifolds 
of a general conservative time-independent flow field F(x). 
Different from the procedure, our robots will solely rely on 
information that can be gathered via local sensing and shared 
through the network. In contrast, a straight implementation 
of the PIM Triple Procedure would require global knowledge 
of the structure of the system dynamics throughout a given 
region given its reliance on computing escape times. We 
describe our approach in the following section. 

B. Controller Synthesis 

Consider a team of three robots and identify them as 
robots {L,C,R}. While the robots may be equipped with 
similar sensing and actuation capabilities, we propose a 
heterogeneous cooperative control strategy. 

Let x(0) = [xl (0), x£(0), xjj(0)] r be the initial conditions 
for the three robots. Assume that x(0) lies on the line 
segment Jq where Jq is a saddle straddle line segment and 
{xl(0),x c (0),x^(0)} constitutes a PIM triple. Similar to the 
PIM Triple Procedure, the objective is to enable the robots 
to maintain a formation such that a valid saddle straddle line 
segment can be maintained between robots L and R. Instead 
of computing the escape times for points on Jq as proposed 
by the PIM Triple Procedure, robot C must remain close 
to Bs using only local measurements of the velocity field 




Fig. 1. Three robots tracking Bg in a conservative vector field. The 
blue dash-dot lines represent the robot trajectories, the green dashed line 
represents the saddle straddle line segment J, and p/. and pn denotes the 
target positions for L and R respectively when executing Up and Ua- 

provided by the rest of the team. As such, we refer to robot 
C as the tracker of the team while robots L and R maintains 
a straddle formation across the boundary at all times. Robots 
L and R may be thought of herding robots, since they control 
and determine the actions of the tracking robot. 

1) Straddling Formation Control: The controller for the 
straddling robots consists of two discrete states: a passive 
control state, Up, and an active control state, Ua- The robots 
initialize in the passive state Up where the objective is to 
follow the flow of the ambient vector field. Therefore, V; = 
for i = L,R. Robots execute Up until they reach the maximum 
allowable separation distance dMax from robot C. When ||x, ■ — 
xc|| > dMax, robot i switches to the active control state, Ua, 
where the objective is to navigate to a point p,- on the current 
projected saddle straddle line segment /, such that ||p,— 
Pe|| = dMin and pc denotes the midpoint of J,. When robots 
execute Ua , V; = ||(p; — x,-) — U;|| and 9j(t) = Ok{t) where 
a, is the angle between the desired, (p; — x,), and current 
heading, u,-, of robot i as shown in Fig. [T] In summary, the 
straddling control strategy for robots L and R is given by 

v = f if d M m < ||x ; — x c || < d Max 

' \ ll(Pi' — x i')~ u i| otherwise 

(3a) 

_ _ f if dMin < ||x; -xcll < d M ax (3h) 
' \ a, otherwise 

We note that while the primary control objective for robots 
L and R is to maintain a straddle formation across Bs, robots 
L and R are also constantly sampling the velocity of the 
local vector field and communicating these measurements 
and their relative positions to robot C. Robot C is then tasked 
to use these measurements to track the position of Bs- 

2) Manifold Tracking Control: Let Ui(t), Uc(f), and &c(t) 
denote the current velocity measurements obtained by robots 
L, C, and R at their respective positions. Let d(-, •) denote the 
Euclidean distance function and assume that d(xc,Bs) < £ 
such that e > is small. Given the straddle line segment J t 
such that xi{k) and Xp(k) are the endpoints of J t , we consider 
an e t < e discretization of J t such that x^ — qi < q2 < . . . < 
qM =Xr. The objective is to use the velocity measurements 
provided by the team to interpolate the vector field at the 



points qi,...,qAf- Since |2} has ff 1 continuity and if Xc is 
e-close to Bs, then the point q# = arg max u(qk) T uc(t) 

k=l,...,M 

should be 5-close to Bs where e < 8 < A and A is a small 
enough positive constant. 

While there are numerous vector field interpolation tech- 
niques available [17]— [19], we employ the inverse distance 
weighting method described in [17]. For a given set of veloc- 
ity measurements u,(f) and corresponding position estimates 
x,(f), the velocity vector at some point qj. is given by 

where Wjj = ||x,-(j) — 1 1 2 . Rather than rely solely on the 
current measurements provided by the three robots, it is 
possible to include the recent history of u,(f) to improve 
the estimate of u(q^), i.e., u,-(f — AT), u,-(f — 2 AT), and so 
on, where Ar is the sampling period and i = {L,C,R}. Thus, 
the control strategy for the tracking robot C is given by 

Vc=||[(qB + M B )-Xc]-Uc|| (4a) 

e c = p c (4b) 

where j3c denotes the difference in the heading of robot C 
and the vector (q^ — up) and b > r is a small number. The 
term bup is included to ensure that the control strategy aims 
for a point in front of robot C rather than behind it. As such, 
the projected saddle straddle line segment J t at each time 
step is given by p 6 = qc + buc with J, orthogonal to Bs at 
qc and ||7 r || chosen to be in the interval [2d m in, 2d Max}- 

IV. ANALYSIS 

In this section, we discuss the theoretical feasibility of the 
proposed saddle straddle control strategy. We begin with the 
following key assumption on the robots' initial positions. 

Assumption 1: Given a team of three robots {L,C,R}, 
assume that d(xc(0),Bs) < £ for a small value of e > 0, 
1 1 xl— xc 1 1 = ||x« — x|| = dMin with dtiin > 2r, and robots L 
and R are on opposite sides of Bs- 

In other words, assume that the robots initialize in a 
valid PIM triple formation and their positions form a saddle 
straddle line segment orthogonal to Bs- Our main result 
concerns the validity of the saddle straddle control strategy. 

Theorem 1: Given a team of 3 robots with kinematics 
given by ([TJ and u,- given by Q, the feedback control strategy 
|3]l and Q maintains a valid saddle straddle line segment in 
the time interval [t ,t +At] if the initial positions of the robots, 
x(f), is a valid PIM triple. 

Proof: To show this, we must show that at time t + At, 
robots L and R remain on opposite sides of Bs- Consider the 
rate of change of the following function 

1 T 

H(x l ,Xr) = -(Xl-Xr) (xl-Xr). 

The above expression is simply one half the square of the 
distance between robots L and R. Let J, denote the saddle 
straddle line segment defined by xp(t) and x^(f) at t and 
let p£ be the intersection of J t and Bs- By construction, if 



we linearize |2]) about the point p#, then the Jacobian of 
|2]i at pa will have one positive eigenvalue. Furthermore, the 
linearized system can be diagonalized such that the direction 
of instability lies along J, [20]. Thus, 4-H > in the time 
interval [t,t + At] when V, = in Q. 

When Vi ^ in |3) for i = L,R, f f H < if the robots 
L and R are moving closer to robot C after reaching the 
maximum allowable separation distance. Recall p m i„(Bs) > 
r, the smallest radius of curvature of B$, and dMin > 2r. 
Furthermore, robot C initializes e-close to the boundary and 
Q steers C towards pc on J, where J, is orthogonal to B$ 
at Xc- This ensures that the rate of the change of the radius 
of curvature of the manifold B$ is small enough such that /, 
intersects with B$ only once. Since dmn > 2r, this ensures 
that even if 4zH < 0, the straddling robots never cross the 
boundary as they move closer to the tracking robot. ■ 

The above theorem guarantees that for any given time 
interval [t,t+At] the team maintains a valid PIM triple 
formation. As such, the iterative application of the proposed 
control strategy leads to the following proposition. 

Proposition 1: Given a team of 3 robots with kinematics 
given by Q and u,- given by Q, the feedback control strategy 
|3]l and |4ji results in an estimate of B$, denoted as Bs, such 
that (Bs,B$) L ^ < W for some W > where (v)l 2 denotes 
the inner product (which provides an measure between 
the Bs and Bs curves). 

From Thm. [T] since the team is able to maintain a valid 
PIM triple formation across Bs for any given time interval 
[f J +At], this ensures that an estimate of Bs in the given time 
interval also exists. Applying this reasoning in a recursive 
fashion, one can show that an estimate of Bs can be obtained 
for any arbitrary time interval. The challenge, however, lies 
in determining the bound on W such that Bs results in a good 
enough approximation since W depends on the sensor and 
actuation noise, the vector interpolation routine, the sampling 
frequency, and the time scales of the flow dynamics. This is 
a direction for future work. 



V. RESULTS 



A. Simulations 



We illustrate the proposed control strategy given by ([3]) 
and Q with some simulation results. Fig. 2(a) shows the 
trajectories of three robots tracking a sinusoidal boundary 



while Fig. 2(b) shows the team tracking a ID star-shaped 
boundary. We note that throughout the entire length of the 
simulation, the team maintains a saddle straddle formation 
across the boundary. 

In both examples, u = — aVq> -Wxif where a,b > 
and <p is an artificial potential function such that <p(x) =0 
for all xeB, and <p(x) < for any x 6 R 2 /B*. The vector 
y/ is a 3 x 1 vector whose entries are given by [0, 0, j{x,y)Y 
where y(x,y) is the curve describing the desired boundary 
[15]. Lastly, the estimated position of the boundary is given 
by the position of the tracking robot, i.e., robot C. In these 
examples, we filtered the boundary position using a simple 
first-order low pass filter. 




(a) 



(b) 



Fig. 2. Trajectories of 3 robots tracking a (a) sinusoidal boundary and a 
(b) star-shaped boundary. The red dashed line is the estimated position of 
the desired boundary shown in solid black. The start positions are shown 
by A and the end positions are shown by the circle-enclosed blue triangles. 





(a) 



(b) 



Fig. 3. Trajectories of the 3 robot team tracking a (a) star shape. The red 
dashed line is the estimated position of the desired boundary shown in solid 
black. The start positions are shown by A and the end positions are shown 
by O- 0>) Snapshot of the multi-robot experiment. 



B. Experiments 

We also implemented the control strategy on our multi- 
robot testbed. The testbed consisted of three mSRV-1 robots 
in a 4.8x5.4 meter workspace. The mSRV-1 are differential- 
drive robots equipped with an embedded processor, color 
camera, and 802.1 1 wireless capability. Localization for each 
robot was provided via a network of overhead cameras. 



Fig. 3(a) shows the trajectories of the robots tracking a star 



shaped boundary shown in black. Fig. 3(b) is a snapshot of 
the experimental run. We refer the interested reader to the 
attached multimedia file for a movie of the full simulation 
and experimental runs. 

VI. EXTENSION TO PERIODIC BOUNDARIES 

In this section, we consider the system of 3 robots with 
kinematics given by ([TJ where u, is determined by the wind- 
driven double-gyre flow model with noise 



-7rAsin(7i:- 



y = nAcos(n 



f(x,t) 



)cos(7i:-) - 
s 

s ax 



lix + i]i{t), 
-Hy + T} 2 (t), 



f(x,t) — £sin(«f + y/)x + (1 — 2esin(o)f + \ff))x. 



(5a) 

(5b) 
(5c) 



When e = 0, the double-gyre flow is time-independent, 
while for e ^0, the gyres undergo a periodic expansion 
and contraction in the x direction. In |5]), A approximately 
determines the amplitude of the velocity vectors, (o/2n gives 



* M 



Fig. 4. Phase portrait of the model given by ((5} with A = 10, fl = 0.005, 
e = o, v = 0, / = 0, and s = 50. 



the oscillation frequency, e determines the amplitude of the 
left-right motion of the separatrix between the gyres, y is the 
phase, jj. determines the dissipation, s scales the dimensions 
of the workspace, and rj,(f) describes a stochastic white 
noise with mean zero and standard deviation a = s/2I, for 
noise intensity /. In this work, 7J,-(?) can be viewed as either 
measurement or environmental noise. Fig. [4] shows the phase 
portrait of the time-independent double-gyre model. 

Fig. [5] shows the use of the control strategy ^ and Q 
to track the Lagrangian coherent structures of the periodic 
double-gyre model with noise. As mentioned in Section [I] 
LCS are extensions of stable and unstable manifolds to non- 
autonomous dynamical systems [21]. We note that while 
the control strategy was based on techniques developed for 
time-independent systems, the method performs surprisingly 
well in tracking LCS for slow time-varying systems in the 
presence of noise. Details regarding LCS computation can be 
found in [5] and we refer the interested reader to the attached 
multimedia file for a movie of the full simulation run. While 
the control strategy was developed for static flows, the movie 
shows the robustness of the strategy for tracking LCS in 
time-varying flows. 

VII. DISCUSSION 

In this paper, we have designed a control strategy that 
allows collaborating robots to track coherent structures and 
manifolds on general static conservative flows. In addition, 
we showed how the strategy can be used to track LCS in 
time-dependent conservative flows with measurement noise. 
The saddle straddle control strategy is based on the com- 
munication of local velocity field measurements obtained 
by each robot. Using the local velocity field information 
provided by the two straddling robots (the herders), one 
robot (the tracker), is able to detect the coherent structures, a 
global structure that delineates the phase space into different 
dynamical regions. Our work is novel in that the robots are 
determining the location of a global structure based solely on 
local information, and as far as we know, the sensing of LCS 
in the ocean has never been performed using autonomous 
vehicles. Moreover, only initial state knowledge of the LCS 
is required locally to get an accurate prediction of the global 
structure. 

While the cooperative control strategy was inspired by 
the PIM Triple Procedure, a procedure that relies on the 



computation of escape times which is a global property of 
the system, the controller itself only relies on information 
provided by each robot's onboard sensors. We also note that 
the realization of the control strategy by the team of robots 
can be achieved without the need for global localization 
information. As such, the strategy is a purely local strategy. 
Furthermore, the cooperative control strategy was derived to 
track the manifolds on a static flow, but performs surprisingly 
well at tracking the LCS in the time-dependent double-gyre 
model in the presence of noise. 

Since realistic quasi-geostrophic ocean models exhibit 
double-gyre flow solutions, our first attempt seems to suggest 
that our methods may be and general enough to be applied to 
more complicated models, including multi-layer PDE ocean 
models. As mentioned in Section |IV] the robustness of the 
control strategy is dependent on numerous parameters in the 
system which includes robots' sensing and communication 
ranges, the bounds on the sensor and actuation noise, the 
vector interpolation technique, the sampling frequency, and 
the relative time scales of the AUV dynamics in relation to 
the surrounding flow dynamics. While our initial results sug- 
gest that our approach may be robust enough to measurement 
noise, a more thorough understanding of the sensitivity of the 
proposed strategy to these various parameters is instrumental 
in extending our approach to more realistic ocean models and 
for field deployment. 

VIII. FUTURE WORK 

In recent years, there has been significant interest in the 
use of AUVs to collect scientific data in the ocean to improve 
our ability to forecast harmful algae blooms and weather 
and climate patterns. One drawback to operating sensors in 
time-dependent and stochastic environments like the ocean 
is that the sensors will tend to escape from their monitoring 
region of interest. As such, the ability to identify and 
track Lagrangian coherent structures (LCS) in these dynamic 
environments is paramount in maintaining appropriate sensor 
coverage in regions of interest. Additionally, since LCS have 
been shown to coincide with optimal trajectories in the ocean 
which minimize the energy and navigation time [3], [4], real- 
time knowledge of these "super-highways" is key in planning 
efficient AUVs paths. 

Of particular interest is the extension of our method to 
more realistic ocean models. Specifically, can we extend 
our current cooperative tracking strategy to a swarm of 
heterogenous mobile and stationary sensors? By increasing 
the team size and incorporating both stationary and mobile 
sensing devices, it is possible to refine our tracking tech- 
nique to reveal the coherent structures at various spatial 
and time scales. One immediate direction for future work 
is to investigate how the proposed strategy scales to larger 
team sizes. Second, underwater environments pose unique 
challenges in terms of wireless communications. In general, 
acoustic transmissions generally have low data rates and 
acoustic wave propagation can be further affected by the 
surrounding fluid dynamics [22]. As such, a second direction 
for future work is to investigate how communication delays 




(e) t=3.0 (f) t=3.2 (g) t=3.6 (h) t=4.0 

Fig. 5. Trajectories of the team of 3 robots tracking the Lagrangian coherent structures of the system described by j3J with A = 10, \l = 0.005, e = 0.1, 
y = 0, / = 0.01, and .s = 50. The trajectories of the straddling robots are shown in black and the estimated LCS is shown in white. 



and missed transmissions impact the overall accuracy of 
the tracking methodology. In this work, we assume an 
initial state knowledge of the LCS is required. This initial 
formation may be difficult to achieve without any prior global 
knowledge of the flows. By considering a team of both 
stationary and mobile sensors, one can potentially obtain an 
initial estimate of a local LCS through the stationary sensing 
network which can then be tracked and further refined by 
the mobile nodes. A third direction for future work is to 
determine how one can strategically place a combination of 
mobile and stationary sensors to provide real-time updates 
on the locations of LCS. 

References 

[1] G. Haller and G. Yuan, "Lagrangian coherent structures and mixing 
in two-dimensional turbulence," Phys. D, vol. 147, pp. 352-370, Dec 
2000. 

[2] S. C. Shadden, F. Lekien, and J. E. Marsden, "Definition and properties 
of lagrangian coherent structures from finite-time lyapunov exponents 
in two-dimensional aperiodicflows," Physica D: Nonlinear Phenom- 
ena, vol. 212, no. 3-4, pp. 271 - 304, 2005. 

[3] T. Inane, S. Shadden, and 1. Marsden, "Optimal trajectory generation 
in ocean flows," in American Control Conference, 2005. Proceedings 
of the 2005, 8-10, 2005, pp. 674 - 679. 

[4] C. Senatore and S. Ross, "Fuel-efficient navigation in complex flows," 
in American Control Conference, 2008, june 2008, pp. 1244 -1248. 

[5] E. Forgoston, L. Billings, R Yecko, and I. B. Schwartz, "Set-based cor- 
ral control in stochastic dynamical systems: Making almost invariant 
sets more invariant," Chaos, vol. 21, no. 013116, 2011. 

[6] C. Hsieh, Z. Jin, D. Marthaler, B. Nguyen, D. Tung, A. Bertozzi, and 
R. Murray, "Experimental validation of an algorithm for cooperative 
boundary tracking," in Proc. of the 2005 American Control Confer- 
ence., 2005, pp. 1078-1083. 

[7] S. Susca, S. Martinez, and F. Bullo, "Monitoring environmental 
boundaries with a robotic sensor network," IEEE Trans, on Control 
Systems Technology, vol. 16, no. 2, pp. 288-296, 2008. 

[8] I. Triandaf and I. B. Schwartz, "A collective motion algorithm for 
tracking time-dependent boundaries," Mathematics and Computers in 
Simulation, vol. 70, pp. 187-202, 2005. 



[9] V. M. Goncalves, L. C. A. Pimenta, C. A. Maia, B. Dutra, and G. A. S. 
Pereira, "Vector fields for robot navigation along time-varying curves 
in n-dimensions," IEEE Trans, on Robotics, vol. 26, no. 4, pp. 647- 
659, 2010. 

[10] F. Zhang, D. M. Fratantoni, D. Paley, J. Lund, and N. E. Leonard, 
"Control of coordinated patterns for ocean sampling," Int. Journal of 
Control, vol. 80, no. 7, pp. 1186-1199, 2007. 

[11] K. M. Lynch, P. Schwartz, I. B. Yang, and R. A. Freeman, "Decen- 
tralized environmental modeling by mobile sensor networks," IEEE 
Trans. Robotics, vol. 24, no. 3, pp. 710-724, 2008. 

[12] W. Wu and F. Zhang, "Cooperative exploration of level surfaces 
of three dimensional scalar fields," Automatica, the IFAC Journall, 
vol. 47, no. 9, pp. 2044-2051, 2011. 

[13] J. Spletzer and R. Fierro, "Optimal positioning strategies for shape 
changes in robot teams," in Proc. of the IEEE Int. Conf. on Robotics 
& Automation, Barcelona, Spain, April 2005, pp. 754-759. 

[14] S. Kalantar and U. R. Zimmer, "Distributed shape control of homo- 
geneous swarms of autonomous underwater vehicles," to appear in 
Autonomous Robots (intl. journal), 2006. 

[15] M. A. Hsieh, S. Loizou, and V. Kumar, "Stabilization of multiple 
robots on stable orbits via local sensing," in Proc. of the Int. Conf. on 
Robotics & Automation (ICRA) 2007, Rome, Italy, April 2007. 

[16] H. E. Nusse and J. A. Yorke, "A procedure for finding numerical 
trajectories on chaotic saddles," Physica D Nonlinear Phenomena, 
vol. 36, pp. 137-156, jun 1989. 

[17] J. C. Agui and J. Jimenez, "On the performance of particle tracking," 
Journal of Fluid Mechanics, vol. 185, pp. 447-4 68, 198 7. [Online]. 
Available: http://dx.doi.org/10.1017/S0022112087003252 

[18] C. Marchioli, V. Armenio, and A. Soldati, "Simple and accurate 
scheme for fluid velocity interpolation for eulerianlagrangian compu- 
tation of dispersed flows in 3d curvilinear grids," Computers & Fluids, 
vol. 36, pp. 1187-1198, 2007. 

[19] E. J. Fuselier and G. B. Wright, "Stability and error estimates for 
vector field interpolation and decomposition on the sphere with rbfs," 
SIAM J. Numer. Anal, vol. 47, pp. 3213-3239, 2009. 

[20] P. Hartman, Ordinary Differential Equations, 2nd ed. Society for 
Industrial & Applied Math, 2002. 

[21] S. Kent, "Lagrangian coherent structures: Generalizing stable and 
unstable manifolds to non-autonomous dynamical systems," Tucson, 
AZ, 2008. 

[22] D. Wang, P. F. Lermusiaux, P. J. Haley, D. Eickstedt, W. G. Leslie, and 
H. Schmidt, "Acoustically focused adaptive sampling and on-board 
routing for marine rapid environmental assessment," Journal of Marine 
Systems, vol. 78, pp. 393-407, 2009. 



